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Abstract 

We present an analytic model of thermal state-to-state rotationally inelastic collisions of polar 
molecules in electric fields. The model is based on the Fraunhofer scattering of matter waves and 
requires Legendre moments characterizing the "shape" of the target in the body-fixed frame as its 
input. The electric field orients the target in the space- fixed frame and thereby effects a striking 
alteration of the dynamical observables: both the phase and amplitude of the oscillations in the 
partial differential cross sections undergo characteristic field-dependent changes that transgress into 
the partial integral cross sections. As the cross sections can be evaluated for a field applied parallel 
or perpendicular to the relative velocity, the model also offers predictions about steric asymmetry. 
We exemplify the field-dependent quantum collision dynamics with the behavior of the Ne-OCS(^S) 
and Ar-NO(^n) systems. A comparison with the close-coupling calculations available for the latter 
system [Chem. Phys. Lett. 313, 491 (1999)] demonstrates the model's ability to qualitatively 
explain the field dependence of all the scattering features observed. 

PACS numbers: 34.10.+X, 34.50.-s, 34.50.Ez 

Keywords: Rotationally inelastic scattering, polar molecules, alignment and orientation, Stark effect, models 
of molecular collisions. 



I. INTRODUCTION 



Collisions of molecules in electric, magnetic, or radiative fields are nearly ubiquitous in 
nature as well as in the laboratory. Molecules colliding in the Earth's atmosphere or in 
interstellar space are commonly subjected to magnetic and radiative fields; in the labora- 
tory, collisions in fields appear with particular prominence in stereo dynamics coherent 
control ^, and molecular trapping and cooling j^. Molecular collisions in fields have been 
the subject of a number of theoretical studies, recently reviewed, e.g., in Ref. ^]. However, 
analytic models of such collisions are scarce, and limited to the collision regime near the 
Wigner limit, see, e.g., ref. [sf. Here we present an analytic model of state-to-state rotation- 
ally inelastic collisions of atoms with polar molecules in electric fields. The model is based 
on the Fraunhofer scattering of matter waves and is applicable to collisions at thermal and 
hyperthermal energies. We develop the model for the collisions of closed-shell atoms with 
(linear) or ^H (symmetric-top equivalent) molecules, and compare it with the close-coupling 



calculations of van Leuken et al. [6|, available for the latter sjptem. 

The field-free Fraunhofer model was developed by Drozdov [7] and generalized by Blair 
in the late 1950s to treat inelastic nuclear scattering. The model provided a much-sought 
explanation of the experimentally observed phase shifts between oscillations in the elastic 
and inelastic differential cross sections for the scattering of protons or a particles by medium- 
sized nuclei, later referred to as the "Blair phase rule." In 1984, the field-free Fraunhofer 
model was adapted by Faubel [9f| to account for rotationally inelastic thermal collisions 
between helium atoms and N2 and CH4 molecules. 

In this paper, we extend the model to include the effects of an electrostatic field on ro- 
tationally inelastic scattering of polar molecules by atoms. Within the model, these effects 
arise due to the orientation of the polar molecules in the space-fixed frame and the concomi- 
tant relaxation of the parity selection rule. Although the model - in both its field-free and 
field-dependent incarnation - is only semiquantitative, it readily explains all the features 
found in the state-to-state differential and integral cross sections and in their dependence 
on the strength and direction of the electrostatic field. These features include the phases of 
the angular oscillations in the differential cross section and their characteristic variation as 
a function of the electric field. 

In Section [ni we prepare the soil by introducing the field-free Fraunhofer model of matter- 
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wave scattering. In Sections IIIII and IIV| we extend the Fraunhofer model to account for 
scattering of polar molecules in electric fields. In Section llllt we work out closed-form 
expressions for the partial and total differential and integral cross sections and the steric 
asymmetry for collisions between closed-shell atoms and polar molecules, and apply 
them to the Ne-OCS(^S, J = ^ J') collision system. Section HVl develops the theory for 
collisions between closed-shell atoms and polar ^11 molecules in electric fields and exemplifies 
the results by treating the Ar-NO (^11, J = | ^ J') collision system. The main conclusions 
of this work are summarized in Section |Vl 



II. THE FRAUNHOFER MODEL OF FIELD-FREE SCATTERING 

We first describe the Fraunhofer model of field-free scattering and discuss its validity. 
The model is based on two approximations: 

(i) The energy sudden approximation, which represents the amplitude 

/.^f(^) = (f|/elW|i) (1) 

for scattering into an angle from an initial, to a final, |f), state in terms of the elastic 
scattering amplitude, fei{i^), at fixed values of the internal coordinates. The energy sudden 
approximation is well justified when the collision energy exceeds the spacing of the internal 
states, Ecou > AEint B, 3- 

(ii) The elastic scattering amplitude fei{^) in Eq. ([1]) is replaced by the amplitude for 
Fraunhofer diffraction by an impenetrable, sharp-edged obstacle as observed at a point of 
radiusvector r from the obstacle, see Fig. [TJ In its simplest form, the Fraunhofer diffraction 
amplitude is given by the integral 

f{^) ^ [ e-'^^^^'^'^dR (2) 



where ip is the polar angle of the radius vector R which traces the shape of the obstacle, 
R = |R|, and k = |k| with k the initial wave vector. Relevant is the shape of the obstacle 
in the space-fixed XY plane, perpendicular to k, itself directed along the space-fixed Z- 
axis, cf. Fig. [H The major approximation made in deriving Eq. ([2]) consists in neglecting 



terms non-linear in r. We note that the notion of a sharp-edged obstacle comes c 



the rigid shell approximation. The latter has been widely used in classical 



10|, |l2|, 



ose to 
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quantum and quasi-quantum treatments of field-free molecular collisions, where 
the collision energy by far exceeds the depth of any potential energy well. 



In optics, Fraunhofer (i.e., far-field) diffraction 15| occurs when the Fresnel number is 
small, 

^ = 4 < 1 (3) 
Here a is the dimension of the obstacle, r = |r| is the distance from the obstacle to the 
observer, and A is the wavelength, cf. Fig. [U Condition ([3]) is well satisfied for nuclear 
scattering at MeV collision energies as well as for molecular collisions at thermal and hyper- 
thermal energies. In the latter case, inequality ([3]) is fulfilled due to the compensation of the 
larger molecular size a by a larger de Broglie wavelength A pertaining to thermal molecular 
velocities. 

We note that the Fraunhofer scattering amplitude, Eq. ([2]), is quite similar to the ampli- 
tude for Born scattering [l^. Either amplitude is a Fourier transform of the target's spatial 
characteristic - either its shape or its potential. Both the Fraunhofer and Born amplitudes 



comprise averages of the phase factor, exp(ikR), over the target's surface or volume [17 1. 

For nearly-circular targets, with a boundary R{(p) = Rq + S{(p) in the XY plane, the 
Fraunhofer integral of Eq. ([2]) can be evaluated and expanded in a power series in the 
deformation S{ip), 

m = /o(^) + 6) + 6') + --- (4) 
with fo{^) the amplitude for scattering by a disk of radius Rq 



.,Ji(Mo'») 



and /i the lowest-order anisotropic amplitude. 



f,{^) = — [ 6{^)e~'^''^°^^'°''^d^ (6) 



2vr^o 

where Ji is a Bessel function of the first kind. Both Eqs. and are applicable at small 
values of ?9 < 30°, i.e., within the validity of the approximation sini) ^ 

A key step required to maintain the analyticity of the Fraunhofer scattering amplitude, 
Eq. ([6]), is to present the shape of the atom- linear molecule potential in terms of a series in 
spherical harmonics, 

R\6\^')=Y,E^,YU0\y^') (7) 
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with S^j, the Legendre moments. The polar and azimuthal angles 6^ and pertain to the 
body-fixed frame, defined, e.g., by the target's principal axes of inertia. However, what 
matters is the target's shape in the space fixed frame, see Fig. [1], which is given by 

/3, 7; e,^) = Y, ^.u%Mh)y.p{o. (8) 

KVp 

where 7) are the Euler angles through which the body-fixed frame is rotated relative 

to the space-fixed frame, {9, (f) are the polar and azimuthal angles in the space-fixed frame, 
and ^pj^lajSj) are the Wigner rotation matrices. Clearly, the term with k, = corresponds 
to a disk of radius Ro, 

"0 - % (9) 



Since of relevance is the shape of the target in the XY plane, we set 6' = | in Eq. (IHl). As a 
result, 

6{^) = R{a,(3,r, -Ro = R{v>) -Ro = Y,^.,%Ml3i)y.p{i.^) (10) 

KVp 

By substituting from Eq. (ITOl) into Eq. and evaluating the integral, we obtain the following 
expression for the first-order scattering amplitude, 

A(a,/3,7;^) = ^5^H«.^;,(a/37)i^.p^|p|(W) (H) 

with Ff^p defined by 

{-l)P27r {'-^)H-tr^^^ foT^ + p even and ^>p 

(12) 

elsewhere 



Fkp= < 



For negative values of p, the factor {—i)'^ is to be replaced by i'^. Finally, by making use of 
Eq. ([T]), we obtain the inelastic scattering amplitude as 

/wf(^) ^ (fl/o + = (f|/i|i) = ^ E ^.u{fmu\i)F.MkRo^) (13) 

KUp 

K+p even 
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III. THE FRAUNHOFER MODEL OF ROTATIONALLY INELASTIC SCATTER- 
ING OF POLAR Is MOLECULES BY CLOSED-SHELL ATOMS IN AN ELECTRO- 
STATIC FIELD 



A. The field-dependent scattering amplitude 

When a polar molecule enters an electrostatic field, its rotational states undergo 
hybridization (coherent linear superposition), induced by the interaction of the molecule's 
body-fixed electric dipole moment, /x, with the electric field, e [3]. Because of the 
cylindrical symmetry about the electric field vector, the permanent- dipole interaction couples 
free-rotor basis states, |J, M), with a fixed value of the good quantum number, M, but a 
range of J's. Thus the hybrid wavefunctions take the general form 

\J,M,u) = J2(^Ui^MM) (14) 
J 

where the expansion coefficients ajj^ depend solely on a dimensionless interaction parameter, 

to = fie/B (15) 

which measures the maximum potential energy, /ie, of the dipole in terms of the molecule's 
rotational constant, B. The symbol J denotes the nominal value of J that pertains to the 
field-free rotational state which adiabatically correlates with the hybrid state, 

\J,M,uj -^0) ^ \J,M) (16) 

and fi = \fj,\, e = \e\. 

In order to account for an arbitrary direction of the electric field with respect to the 
initial wave vector k, we introduce a field-fixed coordinate system X'^Y'^Z'^, whose Z^axis 
is defined by the direction of the electric field vector e. The free-rotor states are thus given 
by spherical harmonics whose arguments are the angles 6"^ and in the field-fixed frame, 

\J,M)=YjM{eK^^) (17) 

Apart from possessing a sui generis energy level pattern, the | J, M, oj) eigenstates have 
an indefinite (mixed) parity and are directional, exhibiting a varying degree of orientation, 
which depends on the values of J, M, and uj. In the oriented states, the body-fixed dipole 
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(and thus the internuclear axis) hbrates about the field direction hke a pendulum, and so the 
hybrid states are referred to as pendular. It is the directionality of the pendular states along 
with their mixed parity that enters the field-dependent Fraunhofer model and distinguishes 
it from the field-free model, which assumes an isotropic distribution of the molecular axes 
and a definite parity. The directional properties of pendular states are exemplified in Figure 
Fig. [2], which shows polar diagrams of the field-free and pendular wave functions at = 5. 

Hence the scattering process in the field comprises the following steps: A molecule in a 
free- rotor state | J, M) enters adiabatically the field where it is transformed into a pendular 
state I J, M, uj). This pendular state may be changed by the collision in the field into another 
pendular state, |J',M',u;). As the molecule leaves the field, the latter pendular state is 
adiabatically transformed into a free-rotor state |J',M'). Thus the net result is, in general, 
a rotationally inelastic collision, |J, M) \ J\M'). 

In order to be able to apply Eq. (fT3l) to collisions in the electrostatic field, we have to 
transform Eq. (fT4l) to the space-fixed frame XY Z. If the electric field vector is specified by 
the Euler angles ((^g, 6'^, 0) in the XY Z frame, the initial and final pendular states take the 
form 

|i) = |J,M,u;) = Y,al,,{uj)Y,&-l^,{ip,,e,,Q)Yj^{e,ip) (18) 
J (. 

(fl = (J',M',u;| = Y,b%,{u)Y,&i:i,,{^,,ee,Q)Y],^.{d,^) (19) 
J' €' 

which is seen to depend solely on the angles and ip (and not the angles 6*** and ip^ pertaining 
to the field- fixed frame). 

On substituting from Eqs. (fT8|) and (fT9!) into Eq. (fT3|) and its integration, we obtain a 
general expression for the Fraunhofer scattering amplitude in the field, 

Kp 

K+p even 



JJ' 



X 5^ajA,(^)6j;;,,(a;)-^|^^C'(J/€J';000)C'(J/«J';MAMM'), (20) 



where AM = M'-M and C{Ji, J2, J3; Mi, M2, M3) are Clebsch-Gordan coeffients j20|. Since 
the atom-linear molecule potential is axially symmetric, only the coefficients contribute 
to the scattering amplitude. 



Eq. (|20I) can be simplified by limiting our considerations to special cases. A first such 
simplification arises when we let the initial free-rotor state to be the ground state, | J, M) = 
|0,0). A second simplification is achieved by restricting the orientation of the electric field 
in the space-fixed frame to a particular geometry. 

(i) For an electric field parallel to the initial wave vector, e fl k, we have 9;, — 0, (/^^ — 0. 
As as result, the Wigner matrix becomes ^Tp^/^M^^i^i^)^ which equals a Kronecker delta, 
S-pAM- Hence only the p = —AM' term yields a nonvanishing contribution to the scattering 
amplitude of Eq. (!20|) , 



K+M' even 

X ^'jo{^)KIA^)J^^C{JkJ'; 000)C{JkJ'; OM'M') (21) 
jj' * 

We see that the angular dependence of the scattering amplitude for the parallel case is 
simple, given by a single Bessel function, J\m'\- 

(ii) For an electric field perpendicular to the initial wave vector, £ _L k, we have 9^ — >■ 
|, ips 0. Hence 

Kp 

K+p even 

X a\{u)h';r,A^)J^^C{JU'- 000)CiJKj'; OM'M') (22) 
JJ' * 

where d'Zp m' ^^^1 Wigner rotation matrices. Since the summation mixes different 

Bessel functions (for a range of p's), the angular dependence of the scattering amplitude in 
the perpendicular case is more involved than in the parallel case. 

We note that, unfortunately, the Fraunhofer model does not distinguish between the 
parallel and antiparallel orientations of the field with respect to the initial wave vector, as 
can be seen by substituting ^!!p_AA/(0' 0) = ^pAAf(-l)''"'' into Eq. ([20]). This defect is 
inherent to the Fraunhofer model, since the diffraction occurs on a two-dimensional obstacle 
in the XY plane, which looks the same from either side of the plane, no matter whether 
£ ft k or £ ti k. 
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B. Results for Ne-OCS(^E, J = ^ J') scattering in an electrostatic field 

We now proceed with the presentation of the colhsional model with a concrete coUision 
system in mind, namely He + OCS(^S, J = ^ J'). The OCS molecule has been widely 
used in experiments with helium nanodroplets |2l|. The electric dipole moment /i = 0.709 D, 
rotational constant B = 0.2039 cm~^, and spectroscopic amicability make the OCS molecule 
a suitable candidate for an experiment to test the field-dependent Fraunhofer model for a 

molecule. 

According to Ref. 22j, the ground-state Ne-OCS potential energy surface has a global 
minimum of a depth of -81.26 cm~^. In order to diminish the effect of this attractive well 
in the collision, we choose a collision energy E'cou = 500 cm~^, which corresponds to a wave 
number k = 21.09 A^^. The "hard shell" of the potential energy surface at this collision 
energy, shown in Fig. [3], we found by a fit to Eq. ([7]) for k < 6. The coefficients E^^q obtained 
from the fit are listed in the Table [B According to Eq. iQ , the Soo coefficient determines 
the hard-sphere radius Ro, which is responsible for elastic scattering. 



1. Differential cross sections 



We start by analyzing the field-free state-to-state differential cross section, which is given 



by 



with 



' {kRo)2 .y{j'+\M'\y.{j'-\M'\y. 

An (J' + \M'\)V.(J'-\M'\)V. 



for J' + \M'\ even 



^J'\M'\ = < 



(23) 



(24) 



otherwise 
see Eq. ( fT3|) and Ref. j9j. We see that the state-to-state differential cross section is pro- 
portional to the square of the Sj/q coefficient, which means that the shape of the repulsive 
potential provides a direct information about the relative probabilities of the field-free tran- 
sitions and vice versa. For the Ne-OCS system, the ^2,0 coefficient dominates the anisotropic 
part of the potential, see Table [H As a result, the corresponding J = — > J' = 2 transition 
is expected to dominate the inelastic cross section. 



Recalling the properties of the Bessel functions 2J], we see that for kRo^ > ^ (which 
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corresponds to > J' degrees for the system under investigation), the differential cross- 
section has the following angular dependence: 



cos^ [kRo^ — for M' even 



sin^ (kRo^ - f ) for M' odd 



(25) 



By averaging over M' and taking into account that ^j'\m'\ vanishes for J' + \M'\ odd, we 
obtain the angular dependence of the differential cross-section for a — > J' transition: 



cos^ [kRo'& — f ) for J' even 



,(?9) ~ < 



(26) 



sin^ [kRo^ - f ) for J' odd 

The "phase shift" of | predicted by Eq. fl26p for the oscillations in the differential cross 
sections corresponding to even and odd field-free transitions is shown in Figs. H] and [5l The 
elastic scattering amplitude, given by Eq. ([5]), has a sin^ [kR^-d — ^) asymptote, and so is 
out of phase with even-J'-transitions. This latter effect, which is known as the "Blair phase 
rule," can be also seen in Figs. [Hand [51 



The state-to-state differential cross sections for scattering in a field parallel (||) and 
perpendicular (±) to k are given by 

Xoii'^^(^) = E^oli'!M'(^) (27) 

M' 

with 

(28) 

and are shown for the Ne-|-OCS collisions at e = 50 and 100 kV/cm in Figs. [Hand [51 The 
figures reveal that an electrostatic field on the order of 10 kV/cm dramatically alters the 
cross-sections. In this subsection we only analyze the field-induced "phase shifts" of the 
oscillations, and defer the discussion of the amplitudes to subsection HUB 21 on the integral 
cross sections, to which the amplitudes are closely related. 

(i) For a parallel field, e || k, the differential cross section, Eq. (I28ll . has the same explicit 
angular dependence as for the field-free case, Eq. (!25ll . However, the field suppresses the 
selection rule ([2^ and so the summation in Eq. ([27]) comprises all M'-states. Therefore, 
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j^,(ll>-L) /Q^ 



'^,(11,^) 



the resulting cross section is a field-dependent mixture of the sine- and cosine-contributions 
given by Eq. fl25l) . The angular dependence of the differential cross sections in Fig. H] can be 
gleaned from Eq. (^T^. The first sum in Eq. (^T^ extends over even k, for even M', and over 
odd K for odd M'. Therefore, the S^o coefficients, Table [H determine not only the relative 
contributions of different J' states, but also of different M' states in Eq. (1271) . Since the H20 
coefficient eclipses the others, transitions to even M' states dominate whenever the field is 
high enough, and the field-free cross-section fl26|) has a cos^ i) asymptote. This can be clearly 
seen in Fig. H) for odd J', there is a field-induced phase shift of the differential cross section, 
which is absent for transitions to even J'. 

(ii) For a perpendicular field, e ± k, several Bessel functions contribute to the scattering 
amplitude. However, since the summation in Eq. fl2^ requires that k, + p he even, it is 
the even Bessel functions which, like in the case of a parallel field, can be expected to 
dominate the cross section. Indeed, the cross sections shown in Figs. [Hand [5] for parallel and 
perpendicular fields are, for J' = 1, 2, 3, similar to one another. However, the J = — J' = 4 
differential cross section in the perpendicular field exhibits an additional phase shift. This 
cross section represents a special case as it is not dominated by the S20 coefficient. The S20 
coefficient fails to dominate the J = ^ J' = A cross section because of the selection rule, 
J' = J; J ±2, that the Clebsch-Gordan coefficients C{J2J', 000) impose on the n = 2 term. 
However, the products of the hybridization coefficients, a'jQ{uj)bj*Q{uj) and ajQ{uj)bj*^2o{^) 
that occur in the term are very small, due to a tiny overlap of the ajQ^u) and bj*,Q{uj) 
hybridization coefficents. Therefore, a superposition of both even and odd Bessel functions 
contributes to the J = ^ J' = A differential cross section. A more detailed discussion of 
the overlaps of the hybridization coefficients follows in the next subsection. 

2. Integral cross sections 

The angular range, 'd < 30°, where the Fraunhofer approximation applies the best, com- 
prises the largest impact-parameter collisions that contribute to the scattering the most, see 
Figs, m and O Therefore, the integral cross section can be obtained to a good approxima- 
tion by integrating the Fraunhofer differential cross section, Eq. (!27j) . over the solid angle 

sin 'dd'ddip, 




(29) 
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The integral cross- sect ions thus obtained for the field parallel and perpendicular to the initial 
wave vector are presented in Fig. [61 One can see that the state-to-state cross section for 
the J = — s> J' = 2 collisions steadily decreases with the interaction parameter tu, whereas 
the other state-to-state cross sections show a non-monotonous dependence. These features 
refiect the dependence of the integral cross sections on the final M' states accessed by the 
inelastic collisions, which is shown in Fig. [7] for the £ || k case. The relative weights of 
different M' states contributing to the cross sections are determined by a combination of 
the coefficients of Eq. (fT3|) and the F^^m' factor of Eq. (fT2!) . Since the S20 coefficient 
looms over the rest and the F^^m' factor is only non- vanishing for k + M' even and n > M', 
collisions leading to M' = 0,2 dominate for all J' values, see Fig. [71 

But how to explain the field dependence of the cross sections for the prevalent J = 
0, M = J' , M' = 0, 2 collisions? The answer comes from the realization that the 
dependence of the integral cross sections on u is entirely determined by the coefficients 
aj^(co') and bj'ij^p in the scattering amplitude - both for e \\ k, Eq. f[2Tl) . and e ± k, Eq. 
(1221) . What the field dependence of the coefficients ajQ^u), 6j,g(co'), and bj,Q{uj) looks like is 
shown in Fig. [HI As noted in subsection IIIIB 11 for J' = 1, 2, 3 the scattering amplitude of 
Eq. f[2Tl) is dominated by the S20 coefficient, which entails the selection rule J' = J; J±2 in 
the summation over J and J'. In the field-free case, this selection rule is only satisfied for 
the J = ^ J' = 2 scattering, which indeed governs the field-free collisions. Once the field 
is applied, the "distributions" of the ajo(co') and bjQ^uj) coefficients undergo a broadening, 
see Fig. [HI (b)-(e). For the J' = 2 channel, such a broadening reduces the overlap of the 
corresponding hybridization coefficients and thus diminishes the J = —>■ J' = 2 cross 
section. In contradistinction, the overlap of the hybridization coefficients for J' = 1,3 
increases with uj, resulting in an increase of J = ^ J' = 1,3 cross sections, see Fig. [61 
At even higher u, the spread of the coefficients is so large that the products ajQ{u)bj,Q{Lj), 
corresponding to the selection rule J' = J; J ± 2, become very small, cf. Fig. [Ht^e). As a 
result, the cross-sections for the J' = 1,3 channels decrease. The field dependence of the 
J = J' = 4 channel is less straightforward, since, as outlined above, its cross section is 
governed by S^o coefficients with k ^2. 

A prominent feature in Fig. [61 is the significant infiuence of the orientation of e with 
respect to k on the cross section for the J = ^ J' = 1 channel. As mentioned above, for 
both £ II k and £ _L k, the partial J = 0; M = — J' = 1; M' = cross section provides the 
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main contribution to the J = ^ J' = 1 transition. However, an inspection of Eqs. (12T|) 
and fl22|) reveals that the integral cross sections for the J = 0; M = ^ J' = 1; M' = 
transitions is always greater for e ± k than for e || k, due to the coefficients (i!!po(f )• 

3. Frontal versus lateral steric asymmetry 
We define the steric asymmetry as 

where the cross sections cr||^_L correspond, respectively, to £ || k and e ± k. The dependence of 
the steric asymmetry on the permanent dipole interaction parameter uj is presented in Fig. [91 
One can see that a particularly pronounced asymmetry obtains for the J = 0^J' = 1,4 
collisions. This can be traced to the field dependence of the corresponding integral cross 
sections. Fig. [6l Within the Fraunhofer model, elastic collisions do not exhibit any steric 
asymmetry. This follows from the isotropy of the elastic scattering amplitude, Eq. (l22l) . 
which depends on the radius Rq only: a sphere looks the same from any direction. 

IV. THE FRAUNHOFER MODEL OF ROTATIONALLY INELASTIC SCATTER- 
ING OF SYMMETRIC-TOP-EQUIVALENT LINEAR POLAR MOLECULES BY 
CLOSED-SHELL ATOMS IN AN ELECTROSTATIC FIELD 

A. The field-dependent scattering amplitude 



Here we consider a symmetric top-equivalent linear polar molecule, such as 0H( Hi), 
NO(^Hi ), colliding with a closed-shell atom. We treat the molecule as a pure Hund's case (a) 
species, characterized by a non-zero projection, Q, of the electronic angular momentum on 
the molecular axis, whose definite-parity rotational wavefunction is given by 



|J,M, \n\,e) = ^ 



|J,M, \n\) + e\J,M,~\n\) 



(31) 



V2[' 

where the symmetry index e distinguishes between the members of a given Q doublet. The 
symmetry index takes the value of +1 for e levels and of —1 for / levels. The parity of the 



wave function is equal to e( 



25| 
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The rotational states of a Hund's case (a) molecule with J > and M > can be oriented 
by coupling the opposite-parity members of an Q doublet via the electric-dipole interaction. 
Such a coupling creates precessing states, in which the body-fixed electric dipole moment 
/X precesses about the field vector. As a result, molecular rotation does not average out 
the dipole moment in first order. A precessing state is a hybrid of the two opposite-parity 
members of an f2-doublet, and can be written as 



|J,M, \n\,w) = a{w)\J,M, |fi|,e = -1) + (3{w)\J, M, |fi|,e = 1), (32) 

with w = fie/ A an interaction parameter which measures the maximum potential energy of 
the electric dipole in terms of the fi-doublet splitting, A, for J = \Q\. For a precessing state 
with w ^ 1, the coefficients 10(^)1 = 1/5(^)1 = 2^^, and the mixing of the states within 
an Q doublet is perfect. A less perfect mixing, |a(w)| 7^ |/?('W^)|; obtains when w < 1. The 
wavefunction, Eq. (|32l) . reduces for a precessing state with a perfect mixing to \J,M,\Q\). 
It is the inherent orientation of the precessing states along with their mixed parity that 
enters the Fraunhofer model for the scattering of Hund's ) molecules in an electric 

field. The directionality of the precessing states is illustrated in Figure [TOl We assume the 
hybridization of J-states for a symmetric-top state to be negligible. 

The symmetric top wavefunction is given by a Wigner rotation matrix 



I J, M, n) = oK 7« = 0) (33) 

In analogy with Eqs. (ITSl) and (fT9|) . we transform the wave function to the space-fixed frame, 

^im{^\ 0\ 0) = J2 ^P/(V^- 0)^5(¥', 0, 0) (34) 

For transitions with \Q\ = the initial and final precessing states can, therefore, be 
written as 



+ [-aiw)+Piw)]^^l\^\iip,e,0)\ (35) 
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+ [-a\w) + P'{w)] (v., e, 0) I (36) 

By substituting from Eqs. fl35|) and fl36l) into Eq. f|T3l) . we finally obtain the scattering 
amplitude for inelastic collisions of symmetric-top molecules in precessing states 



ikRo / 2 J + 1 



47r V 2J 



'+ 1 ^ 



AM{'^e,Oe,0)E^oF^pJ\p\{kRol^) 



Kp 

K+p even 



X C{JkJ'] \n\Q\n\)C{JKj';M^MM'){ [a{w)a\w)+ I3{w)(3\w)] 



'-IT 



+ [a{w)(3'{w) + a{w)(3{w)] 



(37) 



We note that if both the initial and final precessing states are perfectly mixed, the term in 
the curly brackets of Eq. ( 1371) reduces to 2(— 1)'^. The scattering amplitudes for different 
orientations of the electrostatic field s with respect to the initial wave vector k are obtained 
from Eq. (137|) by substituting the appropriate values of the angles: 6*^ = 0; y^e = for e || k, 
and 6'e = |; = for £ ± k. Eq. fl371) implies that the integral cross-sections, cf. Eqs. fl27|) 
- ( 129|) . for J J' transitions are the same in the parallel and perpendicular fields. However, 
the partial integral cross sections for J,M — > J', M' transitions do depend on whether the 
field is parallel or perpendicular to k. 



B. Results for Ar-NO(^ni) scattering in an electrostatic field 

We now consider the excitation of NO(J = = ^ J', \Q\,e/f) by colhsions with 
Ar. under conditions similar to those deBned in Refs. k fe* a hexapole state selector 
selects the e = — !(/) state, Eq. ( 13T]) . which adiabatically evolves into a partially oriented 
state, Eq. ( 132|) . when the collision system enters the electric field. The electric field of 16 
kV/cm, directed parallel to the initial wave vector, e \\ k, creates a precessing state, Eq. (l32ll . 
whose hybridization coefficients are a = 0.832 and f3 = 0.555. Next, a collision with an Ar 
atom excites the NO molecule to a final, field-free state, \ J',M', \Q\,e'). The final, excited 
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state is considered to be exempt from any effects of tlie electric field, as its fi-doubling is 
large and hence w < 1. As a result, P'{w) = 1 or a'{w) = 1 for a final state of e or / parity, 
respectively. For a collision so defined, the scattering amplitude, Eq. (P7|) for e || k, takes 
the form 



ikRo 
4n V 2 J' 



j—-[J\AM\ikRo^) Yl ^^oF.AMCiJKj';MAMM')CiJKj'; |fi|0|fi|) 



X < 



f3{w) 



a[w) 



K+AM even 



AJ 



-1)'^+ (-1 



+ a{w) 



-11 



-1) 



AJ 



> (38) 



\AJ 



where the first or second row of the expression in the curly brackets corresponds to a final 
state of, respectively, e or / parity. The coefficients S^o of the Ar-NO interaction potentia 



27( 1 , are listed in Table [B According to ref . 



extracted from the data of Sumiyoshi et al. 
the Ar-NO potential surface exhibits a global minimum of —115.4 cm~^ and, thus, the 
Fraunhofer model should be valid at collision energies i^cou > 400 cm^^. In ref. 1^, the 
rigid shell QQT model was also used at these energies. 



The state-to-state integral cross sections for spin-conserving collisions = \VL\ 



at 



a collision energy o: 
of Refs. 



442 cm~^ are shown in Fig.[TTl along with the close coupling calculations 
6|] and [26|. The analytic Fraunhofer model provides a simple interpretation of the 
features exhibited by the cross sections. 

First, let us consider the field-free case for an initial / state, i.e., for a{w = 0) = 1 and 

f3{w = 0) = 0: 



MM\{kRo'd) ^ '^^oF^^^MC{JW-M/\MM')C{JU'-\n\{)\n\) 

K+AM even 



■1)'^ + (-1 



\AJ 



lAJ 



(39) 



Eq. fl39l) immediately reveals that if the potential energy surface is governed by terms with k 
even, parity-conserving transitions, / — > /, will dominate for A J even, while parity-changing 
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transitions, / — e, will dominate for A J odd. This propensity can be seen in Fig. [TTJ It was 



28l | by a rather involved analysis of the close-coupling matrix 



explained previously in Ref. 
elements. 

The qualitative features of the scattering in an electric field can also be readily explained 
by the Fraunhofer model. If the field is present and the target molecule oriented, both even 
and odd A J's in the curly brackets of Eq. (l38l) contribute to the scattering for any value of k. 
For a potential energy surface governed by even-K terms, the electric field will enhance the 
parity-conserving transitions for AJ odd, and suppress them for AJ even; parity-breaking 
collisions will prevail for AJ even and subside for AJ odd. 

From Fig. [TT] one could see that for A J > 2 the Fraunhofer model yields an integral cross 
section which is significantly smaller than the one obtained from a close-coupling calculation. 
This is supported by the work of Aoiz et al. 23| , who found that the diffractive contribution 



to the differential cross sections of Ar-NO collisions is much greater for AJ = 2 than for 
AJ = 3. 

As noted in subsection IIV Al the integral cross-sections for J —>■ J' scattering are the 
same for e \\ k and £ ± k. Therefore, the Fraunhofer model does not distinguish between 
the two collisional configurations and yields a zero steric asymmetry as defined by Eq. (|30|) 
for symmetric-top-equivalent molecules. 

Eqs. fl38l) and fl39|) also reveal the angular dependence of the scattering. In particular, by 
making use of the asymptotic forms of the Bessel functions 2^, we see that for a potential 
energy surface governed either by even- or odd-n terms, the differential cross sections for 
parity-conserving and parity-breaking transitions will be out of phase. This is illustrated 
for scattering from |J = |, = |,/) to |J' = |, = |,e//) states in Fig. [12] (full 
curves). We also note that the parity-breaking cross section is much smaller than the 
parity-conserving one, since the Ar-NO potential is dominated by even K-terms, cf. Table [H 
When the field is on, the initial parity is no longer defined. Moreover, both even and odd 
Bessel functions J\p\{kRo{}) contribute to the cross section. As a result, the differential cross 
sections corresponding to the final e and / states (dashed curves) become similar to one 
another. 

We see that the field-free differential cross sections, presented in Fig. [121 are qualitatively 
similar to the results of close-coupling calculations presented in Fig. 4 of ref. {29 1, which 
also show a phase shift between parity-changing and parity conserving cross sections. When 
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the field is turned on, the close coupling calculations also reveal that the cross sections 
corresponding to final e-states exhibit a phase shift and become similar to the cross sections 
for the final /-states, see Fig. 8 of ref. 29(|. 



V. CONCLUSIONS 

We made use of the Fraunhofer model of matter wave scattering to treat rotationally 
inelastic collisions of polar molecules in electric fields. So far, we have worked out the model 
for polar molecules in and ^11 states interacting with closed-shell atoms. In accordance 
with the energy sudden approximation, inherent to the Fraunhofer model, the interaction 
must be dominated by repulsion. This limits the applicability of the model to the thermal 
and hyperthermal collision energy range. However, the model is also inherently quantum 
and, therefore, capable of accounting for interference and other non-classical effects. The 
effect of the electric field enters the model via the directional properties of the molecular 
states and their mixed parity induced by the field. Even a small orientation of the molecule 
is found to cause a large alteration of the scattering observables, such as differential and 
integral cross sections. The strength of the analytic model lies in its ability to separate 
dynamical and geometrical effects and to qualitatively explain the resulting scattering fea- 
tures. These include the angular oscillations in the state-to-state differential cross sections 
or the rotational-state dependent oscillations in the integral cross sections as a function of 
the electric field. In the face of the absence of any other analytic model of collisions in fields, 
the Fraunhofer model is apt at providing a touchstone for understanding such collisions. 
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Figures 




FIG. 1: Schematic of Fraunhofer diffraction by an impenetrable, sliarp-edged obstacle as observed 
at a point of radius vector r{X, Z) from the obstacle. Relevant is the shape of the obstacle in 
the XY plane, perpendicular to the initial wave vector, k, itself directed along the Z-axis of the 
space-fixed system XY Z. The angle ip is the polar angle of the radius vector R which traces the 
shape of the obstacle in the X,Y plane and 'd is the scattering angle. See text. 
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FIG. 2: A comparison of the moduli of the free rotor wavefunctions 1 1 J, M = 0) | , panel (a), with the 
moduli of the pendular wavefunctions ||J, M = 0,u> = 5)|, panel (b). Also shown is the direction 
of the electric field vector, e. 
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270° 

FIG. 3: Equipotential line R^{e'') for the Ne-OCS potential energy surface at a collision energy of 
500 cm^^. We note that the equipotential line for the Ar-NO collision system looks similar and is 
not shown. The Legendre moments for either potential energy surface are listed in Table [H 
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(e) J'=A 




FIG. 4: Differential cross sections for the Ne-OCS ( J = ^ J') collisions in an electrostatic field 
e=50 kV/cm (red dashed line) and 100 kV/cm (blue dotted line), parallel to the relative velocity 
vector. The field-free cross sections are shown by the green solid line. The dashed vertical line 
serves to guide the eye in discerning the angular shifts of the partial cross sections. 
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FIG. 5: Differential cross sections for the Ne-OCS ( J = ^ J') collisions in an electrostatic field 
e=50 kV/cm (red dashed line) and 100 kV/cm (blue dotted line), perpendicular to the relative 
velocity vector. The field-free cross sections are shown by the green solid line. The dashed vertical 
line serves to guide the eye in discerning the angular shifts of the partial cross sections. 
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FIG. 6: Partial integral cross sections for Ne-OCS (J = ^ J') collisions in an electric field 
parallel, panel (a), and perpendicular, panel (b), to the initial wave vector. 
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FIG. 7: Partial integral cross sections for Ne-OCS {J = 0, M = ^ J' , M') collisions in an electric 
field parallel to the initial wave vector. The red solid lines show the M'-averaged cross sections for 
the J = — > J' collisions. 
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FIG. 8: Absolute values of the hybridization coefficients Ojo(w) (black dashed line), 6j,q (red solid 
line) and &j,q (blue solid line) for different values of the interaction parameter oj. The dashed 
vertical lines serves as a guide to the eye. See text. 
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FIG. 9: Steric asymmetry, as defined by Eq. ([301) . for Ne-OCS ( J = — > J') collisions. 




FIG. 10: The moduli of the symmetric-top wavefunction for J = M = $7 = ^. Panel (a) shows 
the field- free wavefunctions, Eq. (I3ip . for e = 1 (blue line) and e = — 1 (red line). Panels (b) 
and (c) show the wavefunctions of the precessing states in the field, Eq. (j32p . for, respectively, an 
incomplete [a = 0.832, /3 = 0.555) and perfect (« = /? = mixing of the Q doublet states. See 
text. 
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FIG. 11: Integral cross sections for the excitation of NO(J = \Q\ = ^,/) in collisions with Ar 
to higher rotational levels of the |0| = ^ manifold. Panels (a) and (b) pertain, respectively, to 
parity-conserving and parity-breaking Ar-NO collisions. The results obtained from the Fraunhofer 
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model are shown by blue curves, those obtained from the close-coupling calculations of Ref. 
by red curves. Dashed lines pertain to field-free scattering, solid lines to scattering in an electric 
field e = 16 kV/cm. 
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FIG. 12: Differential cross sections for the excitation of NO(J = = ^,/) in collisions with 
Ar to the J' = |,r2 = ^ state. Panels (a) and (b) pertain, respectively, to parity-conserving and 
parity-breaking Ar-NO collisions. Green lines pertain to field-free scattering, red lines to scattering 
in an electric field e = 16 kV/cm. 
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Tables 



TABLE I: Hard-shell Lcgcndre moments H^o for the Ne-OCS potential at a collision energy of 500 
cm~^ and for the Ar-NO potential at a collision energy of 442 cm~^. 







(A) 


K 


Ne-OCS 


Ar-NO 





14.7043 


11.0407 


1 


-0.0968 


0.1744 


2 


0.9455 


0.5757 


3 


0.0540 


0.0040 


4 


-0.0384 


-0.0713 


5 


-0.0131 


-0.0013 


6 


0.0012 


0.0106 
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